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ABSTRACT 

We report on radiative hydrodynamic simulations of M dwarf stellar flares and compare the model 
predictions to observations of several flares. The flares were simulated by calculating the hydrodynamic 
response of a model M dwarf atmosphere to a beam of non-thermal electrons. Radiative backwarming 
through numerous soft X-ray, extreme ultraviolet, and ultraviolet transitions are also included. The 
equations of radiative transfer and statistical equilibrium are treated in non-LTE for many transitions 
of hydrogen, helium and the Ca II ion allowing the calculation of detailed line profiles and continuum 
radiation. Two simulations were carried out, with electron beam fluxes corresponding to moderate 
and strong beam heating. In both cases we find the dynamics can be naturally divided into two 
phases: an initial gentle phase in which hydrogen and helium radiate away much of the beam energy, 
and an explosive phase characterized by large hydrodynamic waves. During the initial phase, lower 
chromospheric material is evaporated into higher regions of the atmosphere causing many lines and 
continua to brighten dramatically. The He II 304 line is especially enhanced, becoming the brightest 
line in the flaring spectrum. The hydrogen Balmer lines also become much brighter and show very 
broad line widths, in agreement with observations. We compare our predicted Balmer decrements to 
decrements calculated for several flare observations and find the predictions to be in general agreement 
with the observations. During the explosive phase both condensation and evaporation waves are 
produced. The moderate flare simulation predicts a peak evaporation wave of ~130 km s~^ and a 
condensation wave of ^30 km s~^. The velocity of the condensation wave matches velocities observed 
in several transition region lines. The optical continuum also greatly intensifies, reaching a peak 
increase of 130% (at 6000 A) for the strong flare, but does not match observed white light spectra. 
Subject headings: hydrodynamics — radiative transfer — stars: atmosphere — stars: chromosphere 
— stars: flare — stars: low mass and brown dwarf 



I. INTRODUCTION 

Solar-type flares are common on cool stars with outer convective envelopes. Extensive observations of flares on 
magnetically active M dwarfs (dMe) in the solar neighborhood have revealed many similarities including the presence 
of enhanced and broadened Balmer emission, increases in optical white light and radio emission, and simil ar time scales 
for the impulsive and gradual phases. A detailed comparison of solar and stellar flares is presented in iHaisch et alJ 
([1991). The generally accepted model of solar flares assumes that reconnecting magnetic fields cause the acceleration 
of electrons (and possibly protons). The accelerated electrons travel down magnetic field lines, impacting in the lower 
atmosphere, producing hard X-ray bremsstrahlung and causing increased optical continuum and line emission. As 
the electrons heat the lower atmosphere, it expands and "evaporates" into the transition region and corona, thereby 
increasing t he thernial sof t X-ray flux. The Neup ert effect describes the relationship between h ard and soft X-rays 
on the Sun llNeuD ertlll968l: iMcTiernan et al ]ll999D. and stellar o bservations of the Neupert effect (jHawlev et al.lll995l: 
iGuedel et alJll99(iHHawlev et al.ll2003HMitra-Kraev et alJl2005j) provide evidence that stellar and solar flare processes 
show similar effects. 

Even though solar and stellar flares are likely to be produced by the same basic mechanisms, the different physical 
conditions on M dwarfs compared to the Sun give rise to a number of observational differences. Flares on dMe stars 
occur more frequently and ar e often much more energetic than solar flares. This is likely the result of the stronger 
magnetic fields on dMe stars l|Johns-Krull fc Valentilll996D which produce more magnetic activity than on the Sun. 
For stellar flares, increases in the optical continuum are very common; this is often taken as the defining observational 
signature of a flaring star. However, detectable elevations in the optical continuum are rare in solar flares, since the 
Sun has much brighter background continuum emission, and therefore only the largest flares produce enough white 
light contrast to be seen. 

On the Svm, impulsive hard X-rays are spatially and temporally correlated with increased white light emission 
(jlludson et al. 1992). By extension, the rapid rise in near ultraviolet and blue emission seen in stellar flares is thought 
to correspond with the solar impulsive phase. Therefore, stellar optical flares are likely to be the analog to solar white 
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light flares. We recently developed models of impulsive solar white light flares l) Allred et alJ 120051 hereafter AOS), 
and in this p aper we apply those techniques t o mo del impulsive M dwarf flares. Previous and recent M dwarf flare 
observati ons llRodono et alJll984HDovle e t al."198 8HPhillir)s et al.l|1988t iHawlev fc Pettersenlll99H lEason et alJll992t 
lAbdul-Azi l et alJll995t iJevremovic et all 1998: H awlev et alJl200a) (hereafter, we refer to IHawlev et al.ll200l as H03) 
are used to provide observational constraints on the models. 

The paper is organized as follows. In Section 2, we briefly outline the method by which the radiative hydrodynamic 
equations are solved. We also describe the procedure used to generate the initial atmosphere and the methods for 
incorporating soft X-ray and electron beam heating into the simulations. In Section 3, we discuss the dynamics of our 
simulations and present detailed line and continuum profiles. In Section 4, we compare the M dwarf models to the 
solar models of A05, while in Section 5 we compare predictions of our models to stellar flare observations. In Section 
6 we present our conclusions. 

2. METHOD OF SOLUTION 

The method we have used is described in detail in A05 and lAbbett fc HawlevI l)1999r ). and we refer the reader to 
those papers for a more detailed discussion. The equations of hydrodynamics, statistical equilibrium and radiative 
transfer are solved on a one-dimensional adaptive grid using the RADYN code (Carlsson & Stein 1994, 1995, 1997). 

The preflare atmosphere was generated in a manner similar to AII99. We assume that the plasma moves along 
field lines in a flux tube which we approximate as a one-dimensional cylinder. Our model flux tub e extends from the 
photosphere to the corona and has a length of 10^ cm. The preflare atmosphere of IHawlev fc Fish er (1992) was used 
as the initial starting state for the atmospheric structure. Constant non-radiative quiescent heating was applied to 
grid zones with photospheric column mass (i.e. column mass greater than 3.16 g cm^^) to balance the energy losses in 
the photosphere. The upper boundary was held at 6 x 10^ K. With these boundary conditions and no external sources 
of heating, the atmosphere was allowed to relax to a state of hydrostatic equilibrium. 

We modeled a flare by calculating the radiative hydrodynamic response to a beam of non-thermal electrons injected 
at the top of the loop. The electron beam heats the stellar atmosphere which increases the soft X-ray, extreme 
ultraviolet and ultraviolet radiation (hereafter, we refer to these collectively as XEUV flux); this further heats the 
atmosphere through radiative backwarming. For the dynamic simulations, we change the upper boundary condition 
to a transmitting boundary and remove its temperature constraint. We note that with no external source of coronal 
heating ou r coronal mo del will cool. However, the time scale for coronal cooling is much greater than the impulsive 
time scale l)Cargil]|ll994|) . As the flare evolves, material flows through the upper boundary, emulating the expansion of 
the flux tube. Since our flux tube has a flxed length of 10^ cm, the loop apex and the upper corona move outside the 
region of computation during the course of the flare simulation. Ou r concern in this pa per is modeling optical and UV 
emission which primarily originate from the lower atmosphere. In I Allred et al.l l)2005j) . we found that chromospheric 
heating from the upper corona was negligible when compared to the energy deposited by the electron beam. Thus, 
the "loss" of the upper corona is not likely to affect the chromospheric emission. 

The non-LTE atomic level populations are obtained for a six level with continuum h ydrogen atom, nine leyel wit h 
continuum helium atom and six level with continuum Ca II ion using the technique of iScharmer fc CarlssonI ((1985). 
We calculate numerous bound-bound and bound-free transitions in detail. These are listed in Tables 1 and 2 of A05. 
Complete redistribution is assumed for all transitions, but the effects of partial redistribution are approximated for the 
Lyman transitions by truncating the line proflles at 10 Doppler widths. Other transitions are included as background 
opacity in LTE. The opacities were obtained from the Uppsala package of Gustafsson ( 1973). We include optically-thin 
radiative cooling from bremsstrahlung and coronal metal transitions using emissivities acquired from the ATOMDB 

database ((Smith et al 2001) . The equations are solved on an adaptive grid using 250 grid points in depth, 5 in angle 

and up to 100 in frequency for each transition. 

Numerous observations of hard X-rays during solar flares indicate the presence of downward directed non-thermal 
electrons impacting in the lower solar atmosphere. Unfortunately, current instruments are unable to detect hard X-rays 
produced from flares on stars other than the Sun. However, as mentioned above, observations of the Neupert effect 
provide justification that the heating mechanisms may be the same. Therefore, in these models we employ an electron 
beam as the source o f heating in th e lower stellar atmosphere. In particular, we use an electron beam with energy 
spectrum obtained bv lHolman et al.l ()2003D for the 23 July 2002 X-class solar flare observed by RHESSI. The injected 
electrons are found to have a double power law energy distribution of the form 



El ((<5„ - 2) - (1^) (^„ -8{)\\ (t) (f ) for E, > Es 



where T is the electron energy flux that enters the magnetic loop; Ec is the cutoff energy below which the X-ray 
emission is assumed to be thermal; and Eb is the break energy where the distribution shifts from spectral index 8i to 
We take (5;, (5„, E^ Eb to be constant with values of 3.0, 4.0, 37 keV and 105 keV respectively corresponding to 
the peak of the 23 July 2002 flare. The rate of energy deposited by the electron beam is modeled using the technique 
of Emslic (1978). 

During flares, the density in the transition region and corona becomes elevated increasing the amount of XEUV 
photons originating from these regions. The outward-directed photons we detect as enhanced emission, while 
the downward-directed photons heat the lower atmosphere through increased photoionizations. We have used the 
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ATOMDB database to determine the thermal volume monochromatic emissivity as a function of wavelength and 
temperature. We include line emissivities calculated for approximately 34,000 transitions at 37 temperature points 
ranging from 10^ K to 10^ K. We have grouped the transitions into 14 wavelength bins with a range of 1 - 2500 A. 
Table 3 in A05 lists the range and emissivity-weighted central wavelengths of each bin. The XEUV backwarming rate 
is calculated from Equation 6 of AOS with the assumption that the energy from the photoionized electrons is entirely 
converted to heat. 

3. FLARE SIMULATIONS 

We carried out two flare simulations, corresponding to moderate beam heating with an electron flux, T = 10^*^ ergs 
cm~^ (the FIO flare) and strong beam heating with T = 10^^ ergs cm~^ s~^ (the Fll flare). The FIO flare was 
evolved for 230 s and the Fll flare for 15 s. In both cases the beam energy spectrum is given by the double power law 
as described in Equation ^ 

3.1. FIO Flare Dynamics 

Figure ^ shows the evolution of the temperature and density stratification, the electron beam energy deposition 
rate, the ionization fraction, and the electron density during the FIO flare. The electron beam initially penetrates 
to a depth of 0.26 Mm, and the surrounding temperature quickly rises in response (the first column of Fig. A 
large hydrodynamic wave forms and begins to carry material from the chromosphere into the transition region. By 
0.25 s, the temperature plateaus at ^^9000 K in the region of beam energy deposition as a result of elevated hydrogen 
radiative cooling. For the next several seconds, the flare evolves slowly as much of the beam energy is radiated away 
by hydrogen emission (the second column of Fig. This results in significantly stronger hydrogen emission lines. 
The transition region continues to fill up with plasma evaporated from the chromosphere and becomes extended in 
height (see panels 3 and 4 of Fig. QJ. By 10.0 s most of the hydrogen in the region of beam heating is ionized, and 
the radiative cooling is no longer able to balance the beam heating. The atmosphere rapidly heats to a temperature 
of 50, 000 K, but by 25.0 s it again plateaus (third column of Fig.QJ. This time the equilibrium is a result of helium 
radiative cooling, and a relatively gentle phase begins as hydrogen and helium radiate away most of the beam energy. 
We note that although the plasma in the extended transition region (z=0.8 - 2.2 Mm in panel 4 of Fig.^ appears to 
have "cooled," when compared to the preflare atmosphere, this is not the case. Material from lower, cooler regions is 
brought upward and heated. For example, the plasma located in this region has a temperature of about 60,000 K at 
75.0 s, but originates in the chromosphere from a height of 0.3 Mm and a temperature of 5000 K. In general, when 
comparing the atmospheric structure during the flare evolution, care must be taken to consider that the plasma has 
moved. 

By 98 s, the helium in the region of beam heating has become almost completely ionized and can no longer effectively 
cool the atmosphere (panel 9 of Fig. . The atmosphere responds by explosively heating. The temperature quickly 
rises (panel 2 of Fig. |2J), and a large hydrodynamic wave is created. We refer to this period as the explosive phase. 
Explosive phase dynamics are shown in Figure [2| A high temperature bubble forms in the region near 1.5 Mm (second 
column of Fig. EJ and for the next 130 s the bubble expands as the wave propagates through the atmosphere. This 
stellar flare simulation was not able to track the wave passing through the top of flux tube, and did not attain a 
steady state. Our final atmosphere (last column of Fig. [21) has a denser corona and transition region (panel 8) than the 
preflare atmosphere. We note that much of the original corona has flowed through the original flux tube boundary. In 
A05 we found, for solar flare simulations, that after the wave passed through the boundary of the flux tube the loop 
heated quickly. The solar flare loop top eventually reached a steady state temperature approximately five times the 
preflare value, and we expect a similar result for stellar flare loops. 

The gentle and explosive phase energetics are plotted in Figure El During the gentle phase, the energy deposited 
by the beam is essentially balanced by radiative cooling (panel 1). A small amount of energy goes into increasing the 
helium ionization in the region of 1.0 - 1.5 Mm (see the inset in panel 2). Once the helium is completely ionized, the 
radiative cooling sharply drops, and the explosive phase begins. In the explosive phase, there is a significant portion 
of the beam energy in the region of the expanding high temperature bubble (see panel 2 of Fig. (21) that is not radiated 
away (panel 3 of Fig. ^ . Most of this energy goes into increasing the temperature (panel 4) and expanding the bubble 
(see P d{l/p) in panel 3). 

3.2. Fll Flare Dynamics 

The evolution of the Fll flare is similar to FIO, but proceeds much faster. Figure^lshows the atmosphere at several 
times during this flare. As in FIO, the beam initially penetrates to a height of 0.26 Mm, and in response the atmosphere 
rapidly heats to about 10,000 K. A strong hydrodynamic wave begins to carry material into the transition region. By 
0.03 s, the temperature increase has stabilized as the beam heating is in equilibrium with hydrogen radiative cooling. 
This equilibrium lasts until 0.16 s when most of the hydrogen in the region of beam heating has been ionized. In 
a fashion similar to FIO, the atmosphere again quickly heats but by 0.5 s plateaus at 50,000 K as a result of an 
equilibrium of beam heating with helium radiative cooling (first column of Fig. ^J. By 2.5 s, helium in the region of 
0.3 Mm is completely ionized, and the atmosphere undergoes explosive heating (second column of Fig.^J. As in FIO, 
a high temperature bubble forms and propagates toward the boundary of the fiux tube. A notable difference between 
FIO and Fll is the height at which this happens: ^ 2 Mm in FIO compared to 0.3 Mm in Fll, see panel 2 of Fig. (21 
and panel 2 of Fig. 21 The evolution of Fll is much quicker than FIO, so less material has had time to evaporate into 
the upper atmosphere. The beam therefore, deposits most of its energy lower in the atmosphere in Fll than in FIO. 
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At 15.9 s there is a very steep decrease in temperature at z = 1.5 Mm (panel 4). This is caused by the high level of 
radiative cooling due to increased density at the front of the outward-directed explosive wave. This is likely an artifact 
of the simulation, as strongly blue-shifted transition region material is not seen in flare observations (see S I5.6|I . 

3.3. Line Profiles 

Mass motions during the flare simulations cause significant Doppler shifts and introduce large asymmetries in the 
line profiles. Figure El shows line profiles for four important lines during FIO. The first column (panels 1, 6, 11 and 
16) show that at 1.0 s these emission lines have increased in intensity, but there are no large mass motions since the 
line profiles are relatively symmetric. In the second column, however, the line profiles are significantly distorted as a 
result of large velocities produced by flare heating. The line profiles are especially complicated when they are formed 
in regions of oppositely directed waves. For example, the Ha profile in panel 2 has an enhanced red wing indicating 
downward-directed plasma, but its line center is blue-shifted. The line center is formed higher in the stellar atmosphere 
than the wings, where the plasma is "evaporating" into the transition region. In panel 5, Ha has a very large blue 
wing. This is the result of the high velocity, high density evaporation wave passing through the boundary of the flux 
tube. H/3 responds to flare heating very quickly. By 1.0 s (panel 6), it is already about 150 ti mes b righter than in 
the preflare atmosphere. Also note that H/3 becomes very broad, as discussed in more detail in § 15.31 The He II 304 
line becomes very bright as a result of the extended transition region that forms during the gentle phase (see panel 
4 of Fig. During this phase, it becomes the most intense line in the flaring spectrum. The He II 304 flux reaches 
a maximum at the onset of the explosive phase (Fig.jS) panel 13). Once the explosive phase begins, the atmosphere 
heats beyond 50,000 K (i.e. the temperature of formation of He II) and the intensity decreases (panels 14 and 15). 
In panel 14, a 22 km s~^ wave produced at the onset of the explosive phase causes an enhanced blue peak in the 
He II 304 line proflle. As this wave passes into the higher temperature regions of the corona, it begins to move faster. 
By 230 s, the wave is moving at 129 km s~^ and He II 304 shows a strong blue-shift (panel 15). In contrast, the 
Ca II K line is sensitive to velocities in the chromosphere. In panel 17, Ca ii K is blue-shifted indicating evaporation 
of chromosphcric material into the transition region. The Ca ii K line proflles in panels 18 and 19 are red-shifted 
as a result of chromosphcric condensation waves. At 100 s, a relatively slow (~6 km s~^) downward directed wave 
produced from the initial beam impact passes through the region of Ca ii K formation producing the red shift (panel 
18). At 135 s, a much faster condensation wave (~30 km s~^) produced from the explosive increase in temperature 
passes through this region and Ca ii K is further red-shifted (panel 19). 

3.4. Continuum 

Our flare simulations exhibit strongly enhanced optical continuum emission. Figure El shows optical continuum 
spectra for the preflare atmosphere and a few times during the FIO (panel a) and Fll (panel b) flares. A large Balmer 
jump and noticeable Paschen jump are clearly present in these spectra. The elevated flux peaking at about 6500 A 
is primarily the result of blackbody emission and indicates heating has occurred in the photosphere. The electron 
beam is unable to penetrate deeply enough to directly heat the photosphere; instead the heating in this region is 
primarily due to radiative backwarming from the Balmer and Paschen continua. This is shown in Figure [71 which plots 
the energy deposited in the lower atmosphere as a function of height above the photosphere. Important contributors 
to the gain in internal energy are plotted in Figure [3i showing that energy deposited by radiation is the dominant 
form of flare heating in the photosphere. Figure plots the most important contributors to the radiative heating 
and indicates that Balmer and Paschen continua significantly heat the lower atmosphere. The photosphere reaches 
a maximum temperature increase of 400 K for FIO and 1200 K for Fll. This together with the elevated Paschen 
continuum produces a 32% and a 129% increase in the optical continuum (measured at 6000 A) for FIO and Fll 
respectively. 

4. COMPARISONS TO THE SOLAR FLARE MODELS OF A05 

Compared to the Sun, an active M dwarf has a higher surface gravity, a cooler photosphere, and a hotter corona. 
Figure compares the temperature structure of our M dwarf preflare atmosphere with the solar preflare atmosphere 
used in A05. The higher surface gravity on the M dwarf compresses the atmosphere relative to the Sun, causing the 
chromosphere and transition region to be closer in height to the photosphere than in the solar model. The density of 
the M dwarf corona is approximately an order of magnitude greater than the solar corona. The compressed atmosphere 
results in a much narrower range of beam energy deposition (shown in Figure |S1d). 

Many of the features predicted by the solar models of A05 are similar to those in the M dwarf simulations presented 
here. For example, despite differences in atmospheric structure, both types of simulations predict chromosphcric 
evaporation waves of similar speeds. The reason for this is that the beam energy is deposited in each atmosphere at 
nearly the same column mass. Since the same beam energy is assumed in the solar FIO and M dwarf FIO simulations, 
and the mass above the beam impact location is nearly the same, similar velocities should be expected. Another 
important similarity is that both types of simulations exhibit an initial gentle phase followed by explosive eruption 
in the transition region when the balance between flare heating and radiative cooling is broken. However, there is 
significant difference between the M dwarf and solar models in the onset time of the explosive phase. The explosive 
phase starts at 98 s and 2.5 s for the M dwarf FIO and Fll fiares and at 73 s and 1.0 s for the solar FIO and Fll 
flares. The reason for this is that since the M dwarf model has a higher coronal density, it takes longer for the helium 
to be fully ionized so that the beam heating overcomes the radiative cooling. 
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Another difference is the intensity of the emitted radiation. In the solar model, the helium-emitting plateau in the 
transition region becomes much more extended than in the M dwarf models (compare the first panels of Fig. El and 
Fig. 4 of A05) causing a larger He II 304 enhancement in the solar case. This is seen in Figure O which shows light 
curves for the He II 304 line during the FIO solar and M dwarf simulations. In both cases, the He II 304 flux peaks 
at the onset of the explosive phase. In contrast, the hydrogen Balmcr line emission is stronger in the M dwarf flare 
spectrum. Since the Sun has a hotter photosphere, the continuum level is much higher for the solar Balmer lines. The 
Balmer line profiles for the solar and M dwarf cases are similar, but since the line to continuum ratio is higher in the 
M dwarf case, it has more Balmer line emission. The H/? line profile in Figure [TUl illustrates this effect. 

5. COMPARISON TO M DWARF FLARE OBSERVATIONS 
5.1. Continuum 

Numerous photometric observations of flares on active M dwarfs have shown that the continuum spectr um fits well 
to a blackbody of approximately 9000 K l|Hawlev fc FishedllQQi lEason et allll992l: lAlekseev et al.lll99 4: H03). This 
is puzzling because the electron beam, which is the assumed source of flare heating, is un able to deposit sufficient 
energy in the photosphere to heat it to so high a temperature. One explanation, proposed bv lHawlev fc Fish er (1992), 
suggested that XEUV backwarming may be responsible for heating the photosphere during fiares. To test this, we 
included XEUV backwarming as described in § 2. Figure ITTI plots the energy deposited during the FIO flare due to 
XEUV emission and compares it to beam heating. We find that the XEUV heating is about a hundred times smaller 
than the beam heating and does not significantly contribute to heating the photosphere. As discussed in S I3.4I Balmer 
and Paschen continuum backwarming does heat the photosphere, but even this is insufficient to raise the temperature 
to 9000 K. 

Another explanation is that the observed UBVR enhanced spectrum may not actually be the result of blackbody 
emission, but rather may be due to increased Balmer emission from the chromosphere. This explanation is supported 
by these simulations. Despite a lack of direct photospheric heating, when convolved with the UBVR bandpasses these 
simulations do produce a spectrum with the general shape of a blackbody with a temperature of 8000 - 9000 K. 
Figure ITSl shows an example of the convolved spectrum and blackbody fit. In a fashion similar to H03 we include an 
additional measurement in the UV, centered at 1469 A. The spectrum is similar to observed continuum enhancements 
seen during fiares (cf. Figure 10 of H03), and therefore, it is plausible that the ubiquitous white light continuum may 
not be thermal in nature, but actually a result of Balmer and Paschen emission. However, the increased U band seen 
in our simulations is the result of a strongly e levated Balmer jump which has not been observed in stellar fiare spectra 
ijHawlev fc PettersenllT99lt lEason et alJll992|) . 

5.2. Continuum Dimming 

An i nitial dimming in the stellar conti nuum at the onset of a. w hite light flare has been observed by many au- 
thors l|C ristaldi et al.l 119801 : iGiampaoa et al. 198^ iDovle et alJ 1198 8: Hawlev et al. 199j|). In a manner similar to 
lAbbett fc'^awlevl lll99^and A05, our simulations also exhibit initial continuum dimming. Figure ^] shows light 
curves for Balmer and higher order hydrogen continua. The continuum dimming seen here is a result of non-thermal 
ionization produced by the electron beam and can be understood as follows. The beam raises collisional rates in the 
upper chromosphere causing higher population densities for the excited states of hydrogen. This increases the absorp- 
tion probability for Balmer and higher order hydrogen continuum photons, which previously would have escaped and 
been seen as continuum emission. This results in an initial decrease in the continuum intensity. The time between 
flare onset and the re-brightening of the continuum is controlled by the ratio of recombinations to photoionizations. 
As the flare evolves, the increased electron density in the upper chromosphere elevates the recombination rate and 
begins to reduce the population excess of the excited states of hydrogen. This decreases the photoionization rate, and 
the continuum re-brightens. 

5.3. Line Broadening 

During the im pulsive phase the Balmer lines dramatically increase in width, as observed by numerous authors 
(see for example. iDovl e et al.llT988l iPhillins et al.lll98atfHawlev fc Pettersenlll99lHEason et al.lll992H Abdul-Aziz et all 

^95; H03). Figure H"!] shows line profiles computed during FIO for H/3, and H7 and compares them to line profiles 
observed during th e impulsive peak of Flare 8 of H03. Similarly, Figure ^] shows the H7 line compared to the line 
profile obtained bv lHawlev fc PettersenI l|1991f) . In each case the quiescent line emission has been removed from the 
line profile. To aid in measuring the line width, smoothing has been applied to the H03 observed line profile. To 
include the broadened wings, we measure line widths at one tenth maximum. In both the simulation and the H03 flare 
observation, H/3 increases to a width of about 16 A, and H7 reaches a peak line width of 10 A. The Hawlev & Pcttersep 
((1991) flare was much larger than the H03 flare, and in this flare H7 had a peak width of 15 A compared to 13.1 A 
obtained in our Fll simulation. 

The observed line broadening is likely caused by either Stark broadening from the increase d electron de n sity i n 
the flaring atmosphere , or from turbulent and mass motion flows in the flaring atmosphere. IDovle et alJ l)1988D : 
^ason et al. (1992); Ab dul- Aziz et all l|199,'Tt) found that mass flow velocities of ~100 ~ 300 km s~^ would be required 
to produce the observed Balmer line widths. In both our FIO and Fll models, chromospheric velocities are much less 
than these (H7 has a maximum velocity of ~10 km s~^), and we find that in these simulations the line broadening is 
almost entirely due to the Stark effect. 
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5.4. Balmer Decrement 

The Balmer decrement represents the relative energy emitted in various Balmer lines. It is defined as the ratio of 
the excess flare emission in the Balmer lines relative to a fiducial, often taken to be H/3 or H7. Information about the 
electron density in the Bal mer line forming region of th e atmosphere may be inferred, with steeper decrements implying 
smaller electron densities l|Garcfa- Alvarez et al.ll200^ . A simple explanation for this is that as densities increase, the 
lower order lines become thicker and more energy can escape in the higher order lines, making the decrement shallower. 
Since the decrement is a relative measure, it makes comparing flares of different sizes straightforward. In Table ^ 
we present Balmer decrements calculated relative to H7 for FIO and Fll and compare them to several observations. 
We have also included the Ca 11 K decrement. For the FIO and Fll simulations we employed a six level hydrogen 
atom, and therefore we could not directly calculate Balmer lines higher than II7 in detail. To obtain line fluxes for 
higher order lines we passed a "snapshot" of the atmosphere at the peak of Balmer emission to the static radiative 
transfer code, MULTI (Carlsson 198^. In MULTI we used a 13 level hydrogen atom, allowing us to model the Balmer 
lines through II12. Comparing Fll to FIO shows that for a larger flare, II7 becomes relatively stronger than the other 
Balmer lines. All the observed Balmer decrements more closely match Fll than FIO, as expected si nce these flares were 
all relatively large. The decrements for a flare on AD Leo reported bv lHawlev fc PetterseiJ l|1991[ ) are especially close 
to Fll, indicating that the flare chromospheric structure likely resembles the structure predicted by this simulation. 

5.5. Line to Continuum Ratios 

The continuum carries 79% of the total emission in the FIG flare. The lines are relatively stronger in the Fll flare 
with the continuum carrying only 72% of the energy. Both model predictions are considerably smaller than ratios 
typically observed for M dw arf flares. H 03 observed a continuum emission to total emission ratio of ~92% for five 
moderate flares on AD Leo. iRodono et a l. ( 1984) also observed a ratio of 92% for a flare on AD Leo, and for a very 
large flare iHawlev fc PettersenI (11991') obtained a ratio of 86%. Our over-estimation of the relative line strength is 
likely due to the fact that these simulations do not reproduce the featureless ^9000 K white light continuum observed 
in M dwarf flares. If we add a 9000 K blackbody spectrum to our model spectrum, we do obtain a line to continuum 
ratio of more than 90% as expected. 

Another possible reason for the significant line emission from our simulations is the shape of the beam energy 
spectrum. If we are over-estimating the concentration of beam energy that is deposited in the upper chromosphere, 
the excess heating will produce too much Balmer emission. The majority of the beam energy is deposit ed within 
a rela tively narrow region of the chromosphere and causes large increases in the Balmer emission. iMaiias fc Gonie3 
^^1997^ have shown, using a Fokker-Planck description of the electron beam, that the beam energy deposition is over a 
wider region and higher in the atmosphere than is predicted by the analytic expression employed in these simulations. 
Less concentrated heating will likely produce less dramatic Balmer emission and a line to continuum ratio which more 
closely matches observations. 

5.6. UV Line Velocities 

HOB observed numerous transition region and coronal lines with the HST/STIS instrument during flares on AD 
Leo. They found that many of the lines formed in the transition region were red-shifted during the impulsive phase 
of the flare, and interpreted this as a n example o f chromospheric co ndens ation. To compare our simulations to these 
observations, we used the CHIANTI (|Dere et al.lil997i fYoung et alj 2003) database to produce synthetic line profiles 
for the lines observed by HOB. In our flare simulations, we find multiple regions of the atmosphere with typical 
transition region temperatures (~10^). For example, at 150 s in FIO (panel 3 of Fig. (S)) there are three regions with 
this temperature {z = 0.5, 3.1 and 3.6 Mm). The outer two {z = 3.1 and 3.6 Mm) are traveling upward producing blue- 
shifted emission and the inner one (z = 0.5 Mm) is moving deeper into the atmosphere causing red-shifted emission. 
In Figure ITCl we compare the velocity of the downward directed plasma with the velocities obtained during Flare 8 of 
II03. Both the simula tion and obse rvations exhibit velocities of 20 - 40 km s~^ in agreement with the condensation 
velocities predicted bv lFisherl (|1989D . 

6. CONCLUSIONS 

We have constructed radiative hydrodynamic simulations of a flaring M dwarf atmosphere and used them to explore 
the flare-induced optical and UV emission. We find that the dynamics of the impulsive flare naturally divide into two 
phases, an initial gentle phase followed by a period of explosive increases in temperature and pressure. The explosive 
wave front creates a high temperature bubble which expands until the wave passes through the boundary of the flux 
tube. The atmosphere then attains a steady state with a corona and transition region that are denser than the preflare 
atmosphere. 

Line emission is greatly enhanced in both the moderate (FIO) and strong (Fll) fiare simulations. As chromospheric 
material is evaporated into the transition region, the He II 304 line is especially enhanced, becoming the strongest 
line in the spectrum. Large hydrodynamic waves produced by the flare heating cause significant Doppler shifting of 
the line profiles. The He II 304 line exhibits a peak upward velocity of 129 km s~^, and the Ca 11 K line has a peak 
downward velocity of 30 km s~^. We compared the predicted velocity of our downward-directed condensation wave to 
velocities measured in transition region lines by H03, and found that both models and observations exhibited velocities 
between 20 - 40 km s~^. 

Stark broadening caused by the increased electron density in the flaring atmosphere makes the Balmer lines very wide. 
We compared our FIO Balmer profiles to the observations in H03, and found that for both observations and simulations 
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H/3 had a width of 16 A, and H7 had a width of 10 A. We also found that for the stronger iHawlev fc PetterserJ l)1991f) 
flare, the H7 line was similar in width and shaper to the Fll H7 line profile. 

Our simulations predict elevated optical continuum emission. A large Balmer jump in the flare spectrum radiatively 
backwarms the atmosphere below the temperature minimum region and heats the photosphere by about 400 K in the 
FIG case and 1200 K for Fll. With the elevated Paschen continuum this produces a 32% increase in the continuum 
level for FID and a 129% increase for Fll. The predicted photospheric temperature increase is not large enough to 
explain the observed blackbody continuum of ^^9000 K. Our predicted spectrum, when convolved with the UBVR 
filters, does produce a spectrum which looks like a 9000 K blackbody, but this depends on the presence of a large 
Balmer jump which has not been seen in M dwarf flares. The optical continuum exhibits an initial dimming which 
also has been observed in M dwarf flares. In our simulations, the dimming is a result of non-thermal ionizations due 
to the electron beam causing an over-density of excited states of hydrogen. 

We have successfully modeled velocities, line broadening, and Balmer decrements. While broad band colors are 
matched by our spectrum convolved with filter curves, the continuum does not match in details, predicting a large 
Balmer jump and a noticeable Paschen jump which are not observed. The formation of the strong, featureless white 
light continuum during stellar flares, therefore, remains a mystery. 

Although we have taken great care to make these simulations as accurate as possible, we have found it necessary 
to make several limiting assumptions. The atmosphere is assumed to be one-dimensional, plane-parallel, and aligned 
with the magnetic field. This neglects the possibility that radiation can escape through the sides of the flux tube. 
The ion and electrons temperatures are assumed to be the same. We have also assumed complete redistribution holds 
for all transitions. By neglecting MHD effects we have implicitly assumed that that magnetic pressure is much higher 
than the gas pressure and effectively confines the plasma to the flux tube. 

We plan to make several improvements to these models. As mentioned above the analytic treatment of lEmslid 
l)1978t) results in the beam energy being deposited in a relatively narrow region of the atmosphere causing steep 
electron density gradients and much stronger than observed Balmer emission. We plan to incorporate a model of the 
electron beam using Fokker-Planck kinetic theory. Preliminary analysis shows a wider range of energy deposition with 
a larger portion of the beam deposited in higher regions of the atmosphere than predicted by Emslie (1978). This will 
also provide more direct heating to the corona. With these changes, we hope to be able to model a larger flare loop, as 
is probably more appropriate for the stellar flare case. This may allow the stellar flare explosive wave to pass through 
the top of the loop, and thereby allow the flare loop to obtain a steady state as in the solar case. 

This work has been partially funded by NSF grant AST02-05875 and HST grants AR-10312 and GO-8613. The 
computations presented here were carried out on the Astronomy Condor Network at the University of Washington, 
and we would especially like to thank John Bochanski for contributing many hours of computer time. 
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FIO at t= 0.10 s F10att=3.1s FIO at t= 25.0 s FIO at t= 75.0 s 
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Fig. 1. — The stellar atmosphere at four times during the initial phase of the FIO flare. The top row shows the log of the temperature, 
T (left axis) and velocity (right axis) as a function of height compared with the preflare state (PF). In the middle row the electron density, 
rie (left axis) and beam heating rate, Qe (right axis) are plotted. The bottom row shows the mass density, p (left axis) and hydrogen 
ionization fraction, X, (right axis). Note the change in scale of the horizontal axis in the third and fourth columns. 
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Fig. 2. — The FIO flare atmosphere at four times during the explosive phase. The quantities plotted are identical to Fig. [Tl except that 
the He III fraction is plotted in the bottom row rather the hydrogen ionization fraction. 



10 



AUred, Hawley, Abbett, Carlsson 



-^3 

73 



^ t - — '- Total ^ ^ 
2r_._pd(l/p) 
Radiative 

Conductive 

: Beam 

1 r 



- -1^ 



-1 r 

-2 7 

2 r" 



Time: 75.0 s 



I-total' ' 
-Pd(l/p) 

- - Radiative 

- Conductive 

- Beam 



-1 r 



-2 7 



Time: 110.0 s 




Total 
Thermal 
He ionization 
Other ionization 



(2) 



0.02 
0.01 



Time: 75.0 s 



o.on j . 

1.00 1.25 



1.50 



total 
Thermal 
He ionization 
Other ionization 



(4) 




Time: 110.0 s 



0.5 



1.0 



1.5 2.0 0.5 

Height (Mm) 



1.0 



1.5 



2.0 



Fig. 3. — Energetics for the FIO flare during the gentle phase (panels [1,2]) and the explosive phase (panels [3,4]). The atmospheric 
structure for panels [1,2] is shown in the last column of Figure Q and the structure for panels [3,4] is shown in the second column of 
Figure|^ Important terms in the energy conservation equation are shown in panels [1,3], and contributions to the change in internal energy 
are shown in panels [2,4]. The inset in panel [2] shows a closeup view of the rise in internal energy due to increased He ionization. 
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Fllatt=0.6s Fllatt=2.8s Fllatt=6.0s Fll at t= 15.9 s 
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Fig. 4. — The flare atmosphere at four times during Fll. The quantities plotted are identical to Fig. Ill except that the He III fraction is 
plotted in the bottom row rather the hydrogen ionization fraction. 
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Fig. 5. — Line profiles are plotted for Ha, H/3, He II 304, and Ca ii K at five times during FIO. The smaller line profile in each panel is 
the quiescent line profile. The horizontal dotted line indicates the flare continuum level, and the vertical dotted line is line center. 
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Fig. 6. — Continuum spectra are plotted for the preflare atmosphere and a few times during the FIO (a) and Fll (b) flares. The spectra 
exhibit a large Balmer jump (3646 A), a noticeable Paschen jump (8205 A) and increased blackbody emission. 
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Fig. 7. — Internal energy gain in the lower atmosphere as a function of height above the photosphere, (a) The solid line represents 
the total internal energy gain. The dotted line is the energy deposited by the beam. The dashed line is P dV work, and the dot-dashed 
line is tiic radiative flux. The iicating from radiation dominates in the photosphere, (b) Contributions to the total radiative flux. The 
Balmer, Paschen, Bracket and Pfundt continua are represented by the solid, dotted, dashed, and dot-dot-dot-dashed lines respectively. The 
Balmer and Paschen continua are the dominant sources of heating in this region. Lyman continuum does not significantly contribute to 
the energetics this low in the atmosphere and is not included in the plot. 
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Fig. 8. — (a) Temperature structure of the preflare M dwarf (solid line) and solar (dashed line) atmospheres, (b) Initial beam energy 
deposition rate for the M dwarf (solid line) and solar (dashed line) cases. 
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Fig. 9. — Light curves of the He II 304 Une are plotted for the FIO M dwarf simulation presented in this paper and the FIO solar model 
of A05. The explosive phase begins earlier and there is more total emission in the solar case. 
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Fig. 10. — H/3 line profiles at 100 s for the FIO M dwarf simulation and the FIO model of AOS. The line profiles are similar in shape and 
maximum intensity, but since the level of the continuum is higher in the solar case there is less total line emission. 
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Fig. 11. — Energy deposited during the FIO flare through XEUV backwarming and electron beam heating is plotted as a function of 
atmospheric depth. The XEUV heating accounts for about 1% of the total energy deposited. 
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Fig. 12. — The flare spectrum convolved with the UBVR filters and the corresponding blackbody fit are plotted. The UV point is 
obtained by averaging the spectrum over three 30 A regions centered at 1469 A. The spectrum is best fit with a temperature of 8855 K. 



20 



Allred, Hawley, Abbett, Carlsson 




Fig. 13. — Fluxes for the Balmer, Pasehen, Bracket and Pfundt continua are plotted as a function of time. This plot illustrates the initial 
continuum dimming seen in the higher order hydrogen continua. 
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Fig. 14. — Normalized line profiles for H/3 and H7 obtained from the FIO simulation are compared with observed line profiles at the peak 
of Flare 8 in H03. The line widths at 0.1 maximum are indicated by the dashed lines. The top row shows synthetic profiles from FIO at a 
characteristic time during the initial phase (64 s). The plots on the bottom row show flare spectra (quiescent emission has been removed) 
obtained from the peak of flare 8 from H03. 
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Fig. 15. — Fll H7 line profile compared to the |Hawlev & Pettcrsen (1991) AD Leo flare observation. The simulated line profile is 
taken at the time of peak H7 emission, and the observed line profile is taken at the peak of the flare (542 s). In each case, t he continuum 
emission has been removed and the line proflle has been normalized to the line maximum. The lHawlev Ik Pettersei] 11991]) line width at 
0.1 maximum is indicated by the dotted line. 
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Fig. 16. — FIO plasma velocities are plotted and compared to velocities obtained during Flare 8 of H03. The solid lines plot velocities of 
the plasma at the temperature indicated, and the crosses mark the observed velocities. The labels in each panel indicate the ion measured 
and its formation temperature. The time axis of Flare 8 has been shifted by 7.0 s. 
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TABLE 1 

Balmer Decrement 



Source 


H/3 


H7 


H5 


H8 


H9 


Ca II K 


FIO Flare 


1.28 


1.00 


0.96 


0.82 


0.74 


0.96 


Fll Flare 


1.19 


1.00 


0.89 


0.64 


0.53 


0.12 


AD Leo (HP91)'' 


1.24 


1.00 


0.85 


0.64 


0.47 


0.22 


AD Leo (R84)'' 




1.00 


0.85 


0.53 


0.39 


0.15 


UV Get (PSS)'^ 




1.00 


0.69 


0.47 


0.42 


0.83 


UV Get {E92)<i 




1.00 


0.96 


0.46 


0.34 


0.09 


YZ CMi<= 




1.00 


0.69 


0.54 


0.45 




EZ Aqr* 


1.38 


1.00 


0.83 


0.54 


0.36 




Solar M-class Flare^ 


1.74 


1.00 


0.80 


0.63 







^Impulsive phase. (Hawley fc Pette rsenlll991]) 

'^Impulsive phase. (Rodono et al. 19 8^) 

'^Flare maximum. (Phillips ct al. 19831 

"^Impulsive phase. I Easor^taL 1992) 

•^Flare average. iDo vle et al.lll988 ) 

'Flare maximum o f Flare E. (Jevremovic e t al. 1998) 

^Flare maximum. IJohns-KruU et al.llT997^ 



